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ABSTRACT 


This thesis presents the on-going development of methodology for a probabilistic 
material strength degradation model. The probabilistic model, in the form of a postulated 
randomized multifactor equation, provides for quantification of uncertainty in the lifetime 
material strength of aerospace propulsion system components subjected to a number of diverse 
random effects. This model is embodied in the computer program entitled PROMISS, which 
can include up to eighteen different effects. Presendy, the model includes four effects that 
typically reduce lifetime strength: high temperature, mechanical fatigue, creep and thermal 
fatigue. Statistical analysis was conducted chi experimental Inconel 718 data obtained from the 
open literature. This analysis provided regression parameters for use as the model's empirical 
material constants, thus calibrating the model specifically for Inconel 718. Model calibration 
was carried out for four variables, namely, high temperature, mechanical fatigue, creep and 
thermal fatigue. Methodology to estimate standard deviations of these material constants for 
input into the probabilistic material strength model was developed. Using the current version 
of PROMISS, entided PROMISS93, a sensitivity study for the combined effects of mechanical 
fatigue, creep and thermal fatigue was performed. Results, in the form of cumulative 
distribution functions, illustrated the sensitivity of lifetime strength to any current value of an 
effect In addition, verification studies comparing a combination of mechanical fatigue and 
high temperature effects by model to the combination by experiment were conducted. Thus, 
for Inconel 718, the basic model assumption of independence between effects was evaluated. 
Results from this limited verification study strongly supported this assumption. 
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NOMENCLATURE 


Ai current value of the effect 
Aju ultimate value of the i* effect 
A» reference value of the i** 1 effect 
ai i^ 1 value of the empirical material constant 

b fatigue strength exponent 

c fatigue ductility exponent 

E modulus of elasticity 

K ' cyclic strength coefficient 
n number of effect product terms in the model 
n ' cyclic strain hardening exponent 
N current value of mechanical fatigue cycles 

N ’ current value of thermal fatigue cycles 

Np number of mechanical fatigue cycles to failure 

N'p number of thermal fatigue cycles to failure 
2N’p number of thermal fatigue reversals to failure 
Nu ultimate value of mechanical fatigue cycles 

N'u ultimate value of thermal fatigue cycles 

No reference value of mechanical fatigue cycles 
N’o reference value of thermal fatigue cycles 

q material constant for temperature 

R 2 coefficient of determination 

s material constant for mechanical fatigue cycles 
S current value of material strength 

So reference value of material strength 
T current value of temperature 

Tu ultimate value of temperature 

To reference value of temperature 

t current value of creep time 

4? number of creep hours to failure 
tu ultimate value of creep time 

to reference value of creep time 

u material constant for thermal fatigue cycles 
v material constant for creep time 



NOMENCLATURE (continued) 


Aep/2 plastic strain amplitude 

Aet/ 2 total strain amplitude 

Ao/2 stress amplitude 

c'f fatigue ductility coefficient 

p mean 

a standard deviation 

<j'f fatigue strength coefficient 
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CHAPTER 1 
INTRODUCTION 


Probabilistic methods, for quantifying the uncertainties associated with the design and 
analysis of aerospace propulsion system components, can significantly improve system 
performance and reliability. The reusability and durability of aerospace components are of 
prime interest for economical, as well as, safety related reasons. Life cycle costs including 
initial design costs and field replacement costs of aerospace propulsion system components are 
driving elements for improving life prediction capability. Accurate prediction of expected 
service lifetimes is crucial in the final decision of whether or not to proceed with a particular 
design. Inac curate lifetime strength predictions can result in either a lack of adequate life or an 
overly costly design due to inefficient utilization of material. 

This work is part of a larger effort to develop a probabilistic approach for lifetime 
strength prediction methods [4]. This thesis presents the on-going development of 
methodology that predicts probabilistic lifetime strength of aerospace materials via 
computational simulation. A material strength degradation model, in the form of a randomized 
multifactor equation, is postulated for strength degradation of structural components of 
aerospace propulsion systems subjected to a number of effects. Some of the typical variables 
or effects that propulsion system components are subjected to under normal operating 
conditions include high temperature, fatigue and creep. Methodology to calibrate the model 
using actual experimental materials data together with regression analysis of that data is also 
presented. Material data for the superalloy. Inconel 718, were analyzed using the developed 
methodology. 

Chapters 2 and 3 s ummar ize the theoretical and computational background for the 
research. The above-described randomized multifactor equation is embodied in the computer 
program, PROMISS [6]. This program was developed using the NASA Lewis Research 
Center and the University of Texas System Cray-Y-MP supercomputers. Chapter 4 discusses 
the strength degradation model developed for high temperature, mechanical fatigue, creep and 
thermal fatigue effects, individually. Initial estimates for ultimate and reference values are 
determined using available data for Inconel 7 18. A transformation to improve model sensitivity 
is then discussed. Chapter 5 presents experimental material data for Inconel 718 and displays 
the fate in the form utilized by the multifactor equation embodied in PROMISS. Temperature, 
mechanical fatigue, creep and thermal fatigue data for Inconel 718 are presented. Linear 
regression of the data is performed to provide first estimates of the empirical material constants, 
ai, used to calibrate the model. Additional calibration techniques to improve model accuracy 
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are then discussed. In Chapter 6, methodology for estimating standard deviations of the 
empirical material constants is developed as a means for dealing with limited data. These 
estimated values for the standard deviation, rather than expert opinion, may be used with 
greater confidence in the probabilistic material strength degradation model. Chapter 7 presents 
and discusses cases for analysis that resulted from a sensitivity study for the combined effects 
of mechanical fatigue, creep and thermal fatigue at elevated temperatures. Results, in the form 
of cumulative distribution functions, illustrate the sensitivity of lifetime strength to any current 
value of an effect Chapter 8 presents and discusses model verification studies that were 
conducted to evaluate the ability of the multifactor equation to model two or more effects 
simultaneously. Available data allowed for verification studies comparing a combination of 
mechanical fatigue and temperature effects by model to the combination of these two effects by 
experiment Methodology and results are reiterated and discussed in Chapter 9. Conclusions 
of the current research and recommendations for future research conclude this thesis. The raw 
data for all effects, along with material and heat treatment specifications, are provided in the 
appendix. 
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CHAPTER 2 

THEORETICAL BACKGROUND 


Previously, a general material behavior degradation model for composite materials, 
subjected to a number of diverse effects or variables, was postulated to predict mechanical and 
thermal material properties [7,8,12,13]. The resulting multifactor equation summarizes a 
proposed composite micromechanics theory and has been used to predict material properties for 
a unidirectional fiber-reinforced lamina based on the corresponding properties of the constituent 
materials. 

More recently, the equation has been modified to predict the lifetime strength of a 
single constituent material due to "n" diverse effects or variables [4,5,6]. These effects could 
include variables such as high temperature, creep, mechanical fatigue, thermal fatigue, 
corrosion or even radiation attack. For these variables, strength decreases with an increase in 
the variable [1 1]. The general form of the postulated equation is 


_S_ 

So 



A iU 

A iU 



( 1 ) 


where A if Aju and A;o are the current, ultimate and reference values, respectively, of a 
particular effect; a* is the value of an empirical material constant for the i^ 1 product terms of 
variables in the model; S and So are the current and reference values of material strength. Each 
term has the property that if the current value equals the ultimate value, the lifetime strength will 
be zero. Also, if the current value equals the reference value, the term equals one and strength 
is not affected by that variable. The product form of equation (1) assumes independence 
between the individual effects. This equation may be viewed as a solution to a separable partial 
differential equation in the variables with the further limitation or approximation that a single set 
of separation constants, a;, can adequately model the material properties. 

Calibration of the model is achieved by appropriate curve-fitted least squares linear 
regression of experimental data [18] plotted in the form of equation (1). For example, data for 
just one effect could be plotted on log-log paper. A good fit few the data may be obtained by 
linear regression as shown schematically in Figure 1. Dropping the subscript "i" for a single 
variable, the postulated equation is obtained by noting the linear relation between log S and 
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log [(Au - Aq)/(Au - A)], as follows: 



(2a) 


(2b) 


log S 



Fig. 1 Schematic of Data Illustrating the Effect of One Variable on Strength. 
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This general material strength degradation model given by equation (1), may be used 
to estimate the lifetime strength, S/So, of an aerospace propulsion system component operating 
under the influence of a number of diverse effects or variables. The probabilistic treatment of 
this model includes "randomizing" the deterministic multifactor equation through probabilistic 
analysis by simulation and the generation of probability density function (p.d.f.) estimates for 
life tim e strength, using the non-parametric method of maximum penalized likelihood [19,21]. 
Integration of the probability density function yields the cumulative distribution function 
(c.d.f.) from which probability statements regarding lifetime strength may be made. This 
probabilistic material strength degradation model, therefore, predicts the random lifetime 
strength of an aerospace propulsion component subjected to a number of diverse random 
effects. 

The general probabilistic material strength degradation model, given by equation (1), 
is embodied in the FORTRAN program, PROMISS (Probabilistic Material Strength Simulator) 
[6]. PROMISS calculates the random lifetime strength of an aerospace propulsion component 
subjected to as many as eighteen diverse random effects. Results are presented in the form of 
probability density functions and cumulative distribution functions of lifetime strength, S/So 


CHAPTER 3 

PROMISS COMPUTER PROGRAM 


PROMISS includes a relatively simple "fixed" model as well as a "flexible" model. 
The fixed model postulates a probabilistic multi factor equation that considers the variables 
given in Table 1. The general form of this equation is given by equation (1), wherein there are 
now n = 7 product terms, one for each effect listed below. Note that since this model has 
seven terms, each containing four parameters of the effect (A, Ay, Ao and a), there are a total 
of twenty-eight variables. The flexible model postulates the probabilistic multifactor equation 
that considers up to as many as n = 1 8 effects or variables. These variables may be selected 
to utilize the theory and experimental data currently available for the particular strength 
degradation mechanisms of interest The specific effects included in the flexible model are 
listed in Table 2. To allow for future expansion and customization of the flexible model, six 
"other" effects have been provided. 


Table 1 Variables Available in the "Fixed” Model. 


i* Primitive Primitive 
Variable Variable Type 


1 

Stress due to static load 

2 

Temperature 

3 

Chemical reaction 

4 

Stress due to impact 

5 

Mechanical fatigue 

6 

Thermal fatigue 

7 

Creep 


6 


i 
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Table 2 Variables Available in the "Flexible" Model. 


A. Environmental Effects 

1. Mechanical 

a. Stress 

b. Impact 

c. Other Mechanical Effect 

2. Thermal 

a. Temperature Variation 

b. Thermal Shock 

c. Other Thermal Effect 

3. Other Environmental Effects 

a. Chemical Reaction 

b. Radiation Attack 

c. Other Environmental Effect 

B. Time-Dependent Effects 

1. Mechanical 

a. Creep 

b. Mechanical Fatigue 

c. Other Mech. Time-Dependent Effect 

2. Thermal 

a. Thermal Aging 

b. Thermal Fatigue 

c. Other Thermal Time-Dependent Effect 

3. Other Time-Dependent Effects 

a. Corrosion 

b. Seasonal Attack 

c. Other Time-Dependent Effect 


The considerable scatter of experimental data and the lack of an exact description of 
the underlying physical processes for the combined mechanisms of fatigue, creep, temperature 
variations, and so on, make it natural, if not necessary to consider probabilistic models for a 
strength degradation model. Therefore, the fixed and flexible models corresponding to 
equation (1) are "randomized", and yield the random lifetime material strength due to a number 
of diverse random effects. Note that for the fixed model, equation (1) has the 
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following form: 

S/So = f(Aiu> Aj, Aio, Aju, Aj, Ajo, a,,..., A 7 U, A 7 , A 70 , (3) 

where Ai, Aju and Ajo are the current, ultimate and reference values of the i* of seven effects 
as given in Table 1, and ai is the 1 th empirical material constant In general, this expression can 
be written as, 

S/So = f(X0, i = 1,..., 28, (4) 

where Xj represents the twenty-eight variables in equation (3). Thus, the fixed model is 
"randomized" and assumes all the variables, Xi, i = 1,..., 28, to be random. For the flexible 
model, equation (1) has a form analogous to equations (3) and (4), except that there are as 
many as seventy-two random variables. Applying probabilistic analysis [21] to either of these 
randomized equations yields the distribution of the dependent random variable, lifetime material 
strength, S/So- 

Although a number of methods of probabilistic analysis are available, simulation was 
chosen for PROMISS. Simulation utilizes a theoretical sample generated by numerical 
techniques for each of the random variables [21]. One value from each sample is substituted 
into the functional relationship, equation (3), and one realization of lifetime strength, S/So, is 
calculated. This calculation is repeated for each value in the set of samples, yielding a 
distribution of different values for lifetime strength. 

A probability density function (p.d.f.) is generated from these different values of 
lifetime strength, using a non -parametric method, maximum penalized likelihood. Maximum 
penalized likelihood generates the p.d.f. estimate using the method of maximum likelihood 
together with a penalty function to smooth it [19]. Integration of the generated p.d.f. results in 
the cumulative distribution function (c.d.f.), from which probabilities of lifetime strength can 
be directly noted. 

In summary, PROMISS randomizes the following equation: 




0 ) 


There is a maximum of eighteen possible effects that may be included in the model. For the 
flexible model option, they may be chosen by the user from those in Table 2. For the fixed 
model option, the variables of Table 1 are used. Within the product term for each effect, the 
current, ultimate and reference values, as well as the empirical material constant, may be 
modeled as either deterministic, normal, lognormal, or Wiebull random variables. Simulation 
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is used to generate a set of realizations for lifetime random strength, S/So, from a set of 
realizations for the random variables of each product term. Maximum penalized likelihood is 
used to generate the p.cLf. estimate of lifetime strength, from the set of realizations of lifetime 
strength. Integration of the p.d.f. yields the c.d.f., from which probabilities of lifetime 
strength can be ascertained. PROMISS also provides information on lifetime strength 
statistics, such as the mean, variance, standard deviation and coefficient of variation. 



CHAPTER 4 

STRENGTH DEGRADATION MODEL FOR INCONEL 718 


The probabilistic material strength degradation model, in the form of the multifactor 
equation given by equation (1), when modified for a single effect, results in equation (5) 
below. 



Appropriate values for the ultimate. Ay, and reference quantities, Ao, had to be es tima ted as 
part of the initial calibration of the multifactor equation for Inconel 718. Based on actual 
Inconel 718 data, these values were selected accordingly for each effect. 

4.1 Temperature Model 


Equation (5), when modified for the effect of high temperature only, becomes: 



where Ty is the ultimate or melting temperature of the material. To is a reference or room 
temperature, T is the current temperature of the material, and q is an empirical material constant 
that represents the slope of a straight line fit of the modeled data on log-log paper. A logical 
choice for the ultimate temperature value is the average melting temperature (2369 °F) of 
Inconel 718. Therefore, this value was an initial estimate for the ultimate temperature value, 
Ty. An estimate of 75 °F or room temperature was used for the reference temperature value. 
To. Substitution of these values into equation (6a) above results in equation (6b) below. 
Thus, equation (6b) models the effect of high temperature on the lifetime strength of the 
specified material. Inconel 718. 
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4.2 High-Cycle Mechanical Fatigue Model 


Equation (5), when modified for the effect of mechanical fatigue, becomes: 


_S_ = N u -Np 
S 0 L N u-N 


(7a) 


where Nu is the ultimate number of cycles for which fatigue strength is very small. No is a 
reference number of cycles for which fatigue strength is very large, N is the current number of 
cycles the material has undergone, and s is the empirical material constant for the high-cycle 
mechanical fatigue effect. An initial estimate of lxlO 10 was used for the ultimate number of 
cycles, Nu, since mechanical fatigue data beyond this value was not found for Inconel 718. An 
initial estimate of 0.5 or half a cycle was used for the reference number of cycles. No. 
Substitution of these values into equation (7a) results in the high-cycle mechanical fatigue 
model for Inconel 718, as given below by equation (7b). 


s _r io io -o.5 i 5 

S 0 ~ 10 10 — N 


(7b) 


Since the high-cycle fatigue domain is associated with lower loads and longer lives, or high 
numbers of cycles to failure (greater than 10 4 or 10 5 cycles), data consisting of cycle values 
less than 5X10 4 fall into the low-cycle fatigue regime and therefore, may be modeled by a low- 
cycle mechanical fatigue model rather than the high-cycle one presented here. 


4.3 Creep Model 


Equation (5), when modified for the effect of creep, becomes: 



where tu is the ultimate number of creep hours for which rupture strength is very small, to is a 
reference number of creep hours for which rupture strength is very large, t is the current 
number of creep hours, and v is the empirical material constant for the effect of creep. An 
initial estimate of lxlO 6 was used for the ultimate number of creep hours, tu, due to the fact 
that creep rupture life data beyond this value was not found for Inconel 718. An initial estimate 
of 0.25 hours or fifteen minutes was used for the reference number of creep hours, to- 
Substitution of these values into equation (8a) results in equation (8b) below. 


_S_ 

So 


10 6 -0.25 

10 6 -t 


I v 


(8b) 
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4.4 Thermal Fatigue Model 

The fourth and final effect for which Inconel 718 data was obtained is thermal 
fatigue. Thermal fatigue has been extensively discussed in the literature [9, 16, 23]. When 
modified for the effect of thermal fatigue, equation (5) becomes: 



where Nu is the ultimate number of thermal cycles for which thermal fatigue strength is very 
small, N’o is a reference number of thermal cycles for which thermal fatigue strength is very 
large, N' is the current number of thermal cycles the material has undergone, and u is an 
empirical material constant that represents the slope of a straight line fit of the modeled data on 
log-log paper. 

Thermal fatigue is in the regime of low-cycle fatigue (less than 10 4 or 10 5 cycles), 
therefore, an intermediate value of 5x10* cycles was an initial estimate for the ultimate number 
of thermal fatigue cycles, N’u- An initial estimate of 0.5 or half a cycle was used for the 
reference number of cycles, N'o- Substitution of these values into equation (9a) results in the 
thermal fatigue model for Inconel 7 1 8, as given by equation (9b) below. 



4.5 Model Transformation 


In the case of mechanical fatigue, creep and thermal fatigue, the current value and the 
reference value are small compared to the ultimate value. Therefore, regardless of the current 


value used, the term 


Ap ~ A 
Ait An 


remains approximately constant at a value of one. In order to 


sensitize the model for these three effects, the logio of each value was used. As seen in 
Tables 3 through 5, this transformation significantly increases the sensitivity of a product term 
to the data used within it In addition, this transformation results in better statistical linear 
regression fits of the data, as seen later in Figures 6, 9 and 17 of Chapter 5. Hence, the 

general term [ — — was modified to the sensitized form, ^ — ^°S(A) f or 

[Au-AoJ [ log(A u ) - log( A 0 ) J 
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these three effects. The program, PROMISS93, modifies the program, PROMISS, to allow 
for the sensitized form of these three effects. 


Table 3 Non-sensitized and Sensitized Terms for High-Cycle Mechanical Fatigue Data. 

Test Temperature, 

Cycles, 

r (io io )-(n)' 




°F 

N 





75 

105 

0.99999 

0.485388 


10 6 

0.9999 

0.388311 


10 7 

0.999 

0.291233 


10* 

0.99 

0.194155 

1000 

105 

0.99999 

0.485388 


10 6 

0.9999 

0.388311 


10 7 

0.999 

0.291233 


10 8 

0.99 

0.194155 


1200 


10 5 

0.99999 

0.485388 

10 6 

0.9999 

0.388311 

10 7 

0.999 

0.291233 

10 8 

0.99 

0.194155 


Table 4 Non-sensitized and Sensitized Terms for Thermal Fatigue Data. 


Cycles, 

N' 

r(5*10 J )-(N)l 


log(5xl0 4 )-log(N ) 


log^5 x 10 4 j - log(0.5) 

45 

0.999110 

0.609151 

140 

0.997210 

0.510568 

750 

0.985010 

0.364782 

9750 

0.805008 

0.141993 
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Table 5 Non-sensidzed and Sensitized Terms for Creep Rupture Data. 


Test Temperature, Rupture Life, t. 

' K)-(t) ' 


log(l0 6 )-log(t) 

°F Hrs 

(lO 6 ) — (0.25) 


log(l0 6 ) - log(0.25) 


1000 

27.8 

133.2 

256.0 

814.9 

1731.0 

8473.0 
21523.6 

0.99997 

0.99987 

0.99974 

0.99919 

0.99827 

0.99153 

0.97848 

0.69008 

0.58701 

0.54404 

0.46787 

0.41831 

0.31384 

0.25251 

1100 

28.2 

0.99997 

0.68914 


62.0 

0.99994 

0.63732 


151.9 

0.99985 

0.57837 


367.5 

0.99963 

0.52025 


2327.6 

0.99767 

0.39883 


10606.2 

0.98939 

0.29906 


33990.7 

0.96601 

0.22245 

1200 

10.6 

0.99999 

0.75351 


30.8 

0.99997 

0.68334 


150.0 

0.99985 

0.57920 


747.2 

0.99925 

0.47357 


3131.5 

0.99687 

0.37931 


7263.0 

0.99274 

0.32397 


10232.0 

0.98977 

0.30143 

1300 

18.0 

0.99998 

0.71867 


70.5 

0.99993 

0.62887 


182.7 

0.99982 

0.56623 


476.8 

0.99952 

0.50313 


808.0 

0.99919 

0.46843 


2870.7 

0.99713 

0.38503 


6048.0 

0.99395 

0.33601 


CHAPTER 5 

EXPERIMENTAL MATERIAL DATA 


In order to calibrate or anchor the empirical material constants, ai, in the multifactor 
equation to particular aerospace materials of interest, it is necessary to collect experimental data. 
Since actual experiments were not conducted as part of this research project, data for several 
effects were collected from the open literature. 

5.1 Literature Search 

Initially, a computerized literature search of nickel-base superalloys was conducted to 
obtain existing experimental data on various material properties. Useful data on high 
temperature, mechanical fatigue and creep properties were found for several nickel-base 
superalloys [2, 10, 14, 22]. Based on this data, a second computerized literature search of the 
superalloy, Inconel 718, was later performed in an attempt to find additional data, especially 
data on thermal fatigue effects. Efforts were concentrated on this particular superalloy for two 
primary reasons. First, Inconel 7 1 8 was selected as (he initial material to be analyzed due to its 
extensive u tilizati on by the aircraft and aerospace industries owing to its high performance 
properties. Secondly, data on Inconel 718 was far more abundant than for any other 
superalloy. As a result, data for three effects, namely, high temperature, mechanical fatigue 
and creep were readily obtained. Data on thermal fatigue properties, however, was much 
harder to obtain. Therefore, a third computerized literature search for Inconel 718 thermal 
fatigue data was required. This search yielded limited thermal fatigue data for Inconel 7 18. 

5.2 Inconel 718 

Inconel 718 is a precipitation-hardenable nickel-chromium alloy containing significant 
amounts of iron, niobium and molybdenum along with lesser amounts of aluminum and 
titanium. It combines corrosion resistance and high strength with outstanding weldability. 
Inconel 718 has excellent creep-rupture strength and a high fatigue endurance limit up to 
1300 °F (700 °Q. It requires a somewhat complex heat treatment (solution anneal, cool and 
duplex age) to produce its high strength properties. Standard production forms are round, 
flats, extruded section, pipe, tube, forging stock, plate, sheet, strip and wire. Inconel 718 
material in various forms is used in gas turbines, rocket engines (including the space shuttle 
main engine), spacecraft structural components, nuclear reactors, pumps and tooling. In gas 
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turbine engines, for example, components operate under rigorous conditions of stress and 
temperature. The high performance superalloy. Inconel 718, is capable of meeting such 
extreme material requirements. 

5.3 Temperature Data 

The data on high temperature tensile strength properties of Inconel 718 resulted from 
tests conducted on hot-rolled round specimens annealed at 1950 °F and aged. [14]. This data, 
as well as the data on mechanical fatigue, creep, and thermal fatigue strength properties, were 
plotted in various forms, one of which was the same as that used by the multifactor equation in 
PROMISS. The data plotted in Figures 2 and 3 show the effect of temperature on yield 
strength for Inconel 718. Figure 2 displays the raw data, while Figure 3 shows the data in the 
form given by equation (6b). As expected, the yield strength of the material decreases as the 
temperature increases. Linear regression of the data, as seen in Figure 3, produced a first 
estimate of the empirical material constant, q, for the temperature effect This estimated value 
of the material constant, q, is given by the slope of the linear regression fit. As seen by 
Figure 3 and corroborated by the high R 2 (coefficient of determination [3] ) value, this 
temperature data, when modeled by equation (6b), does indeed indicate a good linear relation 
between yield strength and temperature. 
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Fig. 2 Effect of Temperature (°F) on Yield Strength for Inconel 7 1 8. 
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5.4 High-Cycle Mechanical Fatigue Data 

The data on mechanical fatigue strength properties resulted from high-cycle fatigue 
tests conducted on hot-rolled bar specimens annealed at 1750 °F and aged [14]. This data was 
plotted in various forms, including non-sensitized and sensitized model forms. Figure 4 
presents the raw mechanical fatigue data and displays the effect of mechanical fatigue cycles on 
fatigue strength for given test temperatures. As expected, the fatigue strength of Inconel 7 1 8 
decreases as the number of cycles increases. Figures 5 and 6 show the data in the non- 
sensitized form of equation (7b) and the sensitized model form, respectively. Linear regression 
of the data produced first estimates of the empirical material constant, s, for the mechanical 
fatigue effect, as given by the slopes of the linear regression fits. As seen by these regression 
fits in Figures 5 and 6, the R 2 (goodness of fit) values are significantly higher for the sensitized 
model form. 

In reference to Figure 6, the R 2 value corresponding to a temperature of 75 °F is 
significantly lower than the fits calculated at temperatures of 1000 °F and 1200 °F. In addition, 
whereas the slope corresponding to a temperature of 1000 °F is lower than that corresponding 
to 1200 °F, the slope obtained at a temperature of 75 °F (s = 0.37848) is higher than that at both 
1000 °F (s — 0.22348) and 1200 °F (s = 0.35425). This is due to the fact that at certain current 
cycle values, N, the fatigue strength at a temperature of 75 °F is lower than that at 1000 °F. 
Since this phenomenon is highly improbable, the validity of the mechanical fatigue data 
obtained at a test temperature of 75 °F is questionable. Thus, the corresponding mechanical 
fatigue material constant (s = 0.37848) is also questionable. 



Fig. 4 Effect of Mechanical Fatigue (Cycles) on Fatigue Strength for Inconel 718. 
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0.000 0.001 0.002 0.003 0.004 0.005 

LOG[(10 A 10-0.5)/(10M0-N)] 

Fig. 5 Effect of Mechanical Fatigue (Cycles) on Fatigue Strength for Inconel 718. 

(Non-sensitized Model Form) 



0.3 0.4 0.5 0.6 0.7 0.8 


LOG[(LOG (1 0 A 1 0)-LOG(0.5))/(LOG(1 0 A 1 0)-LOG(N))] 

Fig. 6 Effect of Mechanical Fatigue (Cycles) on Fatigue Strength for Inconel 718. 

(Sensitized Model Form) 
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5.5 Creep Rupture Data 

The data on creep rupture strength properties resulted from tests conducted on stress 
rupture test bars annealed at 1800 °F and aged [2]. As with the mechanical fatigue data, this 
data was plotted in various forms. Figure 7 presents the raw creep rupture strength and 
shows the effect of creep time on rupture strength for given test temperatures. Once again, the 
strength of the material decreases as the variable, in this case time, increases. In addition, for a 
given time, t, the rupture strength decreases as the test temperature increases. This 
phenomenon is clearly seen in Figure 7, as well as, by the changing slopes of the linear 
regression fits in Figures 8 and 9. Figures 8 and 9 show the creep data in the non-sensitized 
form of equation (8b) and the sensitized model form, respectively. Linear regression of the 
data produced first estimates of the empirical material constant, v, for the creep effect, as given 
by the slopes of the linear regression fits. As seen by these regression fits in Figures 8 and 9, 
the R 2 (goodness of fit) value is significantly higher for the sensitized model form. 



RUPTURE LIFE, t p (HOURS) 


Fig. 7 Effect of Creep Time (Hours) on Rupture Strength for Inconel 718. 

(Linear Plot) 
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Fig. 8 Effect of Creep Time (Hours) on Rupture Strength for Inconel 718. 
(Non-sensidzed Model Form) 



R A 2 x 0.994 
R A 2 x 0.992 

R A 2 x 0.991 
R A 2 x 0.998 


Fig. 9 Effect of Creep Time (Hours) on Rupture Strength for Inconel 718. 

(Sensitized Model Form) 
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5.6 Thermal Fatigue Data 

Low cycle fatigue produces cumulative material damage and ultimate failure in a 
component by the cyclic application of strains that extend into the plastic range. Failure 
typically occurs under 10 4 or 10 5 cycles. Low cycle fatigue is often produced mechanically 
under isothermal conditions. However, machine components may also be subjected to low- 
cycle fatigue due to a cyclic thermal field. These cyclic temperature changes produce ther mal 
expansions and contractions that, if constrained, produce cyclic stresses and strains. These 
thermally induced stresses and strains result in fatigue failure in the s am e manner as those 
produced mechanically. 

The general model for the thermal fatigue effect uses stress-life (o-N) data obtained 
from experimental strain-life (e-N) data. The thermal fatigue data presented in Table 6 resulted 
from thermomechanical fatigue tests conducted on test bars annealed at 1800 °F and aged [16]. 
The temperature and strain were computer-controlled by the same triang ular waveform with in- 
phase cycling at a frequency of 0.0056 Hz.. The temperature was cycled between a mi nimnrn 
temperature of 600 °F and a maximum temperature of 1200 °F, with a mean temperature of 
approximately 900 °F. This total strain amplitude data and plastic strain amplitude data were 
used to construct the strain-life curves presented in Figure 10. 


Table 6 Thermal Fatigue Data for Inconel 718. 


Cycles to 
Failure 
N'f 

Total Strain 
Amplitude, 

Ae-i/2 

Plastic Strain 
Amplitude, 
A£p/2 

Stress 
Amplitude, 
AO/2 (psi) 

45 

0.0100 

0.0050 

126,500 

140 

0.0075 

0.0029 

116,380 

750 

0.0050 

0.0011 

98,670 

9750 

0.0040 

0.0003 

93,610 


i 
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Fig. 10 Strain-life Curve for Inconel 718. 


By equation (10), the stress amplitude, Ao/2, was calculated using total and plastic 
strain amplitudes, Ae T /2 and Aep/2, respectively, along with an average value of E=25xl0 6 psi 
(modulus of elasticity for Inconel 718 at 900 °F [14J ). 



The resulting stress amplitude data were then plotted against the plastic strain amplitude data to 
produce the cyclic stress-strain curve shown below in Figure 11. 



Fig. 1 1 Cyclic Stress-Strain Curve for Inconel 7 1 8. 



24 


Using power law regression techniques [1] and the data in Table 6, thermal fatigue 
properties for Inconel 718 were calculated. These properties were calculated and compared 
with known established values in order to check the validity of the data. The plastic portion of 
the strain-life curve (Figure 10) may be represented by the following power law function: 

-^ £ - = e F (2N F ) , ( 11 ) 

where Agp/2 is the plastic strain amplitude and 2N' F are the reversals to failure. A power law 
regression analysis of the data yielded two thermal fatigue properties, namely, the fatigue 
ductility coefficient, £ F , and the fatigue ductility exponent, c. These two properties are 
indicated graphically, along with their coefficient of determination, R 2 , in Figure 12. 
Regression statistics, such as R 2 , were obtained to indicate whether or not a power law 
representation of the relationship between plastic strain amplitude and reversals to failure was 
appropriate. As confirmed by the high R 2 value in Figure 12, the power law function of 
equation (1 1) well represents the relationship between A£p/2 and 2N’ F . 



Fig. 12 Regression of Equation (1 1) Data Yielding Fatigue Ductility Coefficient, e’p, 

and Fatigue Ductility Exponent, c. 
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The following power law function was satisfactory for expressing the cyclic stress- 
strain relationship of the data presented in Figure 1 1 : 



Regression analysis of this data yielded two more thermal fatigue properties, K , the cyclic 
strength coefficient and n', the cyclic strain hardening exponent. These two properties are 
indicated graphically, along with their coefficient of determination, R 2 , in Figure 13. 



Fig. 13 Regression of Equation (12) Data Yielding Cyclic Strength Coefficient, K', 
and Cyclic Strain Hardening Exponent, n\ 

The following power law function was used to approximate the relationship between 
stress amplitude and reversals to failure: 

^ = o' p ( 2N f )\ (13) 

Regression analysis of this data yielded two more thermal fatigue properties, o’ f, the fatigue 
strength coefficient and b, the fatigue strength exponent. These two properties are indicated 
graphically, along with their coefficient of determination, R 2 , in Figure 14. They complete the 



LOG STRESS AMPLITUDE (PSI) 


set of thermal fatigue material properties calculated. The complete set of properties are given in 
Table 7, along with accepted ranges for the exponents [1]. 



Fig. 14 Regression of Equation (13) Yielding Fatigue Strength Coefficient, op, 

and Fatigue Strength Exponent, b. 


Table 7 Thermal Fatigue Material Properties for Inconel 718. 


Material Property 

Calculated 

Value 

Accepted 

Range 

Fatigue Ductility Coefficient, e'p 

-1.2637 

(0.0545) 


Fatigue Ductility Exponent, c 

-0.5279 

-0.5 to -0.7 

Cyclic Strength Coefficient, K’ 

5.3416 

(219,584 psi) 


Cyclic Strain Hardening 
Exponent, n’ 

0.1089 

0.10 to 0.25 

Fatigue Strength Coefficient, o’p 

5.2031 

(159,625 psi) 


Fatigue Strength Exponent, b 

-0.0572 

-0.05 to -0.07 
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The thermal fatigue stress- life (a-N) data were plotted in various forms. Figure 15 
presents the thermal fatigue data and displays the effect of thermal fatigue cycles on stress 
amplitude at failure (i.e., thermal fatigue strength) for a mean thermal cycling temperature of 
900 °F. As expected, the thermal fatigue strength decreases as the number of cycles increases. 
Once again, the data was plotted in both non-sensitized and sensitized model forms to illustrate 
how the sensitized model results in a significant increase in the R 2 (goodness of fit) value. 
Figure 16 presents the data in the non-sensitized form of equation (9b), while Figure 17 shows 
the data in the sensitized model form. Linear regression of the data, as seen in Figure 17, 
produced a first estimate of the empirical material constant, u, for the thermal fatigue effect, as 
given by the slope of the linear regression fit 



Fig. 15 Effect of Thermal Fatigue (Cycles) on Thermal Fatigue Strength 
(i.e. , Stress Amplitude at Failure) for Inconel 718. 



LOQ THERMAL FATIGUE STRENGTH (PSI) LOG THERMAL FATIGUE STRENGTH (PSI) 
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Fig. 16 Effect of Thermal Fatigue (Cycles) on Thermal Fatigue Strength. 
(Non-sensitized Model Form) 



LOG[(LOG(5x1 0 A 4)-LOG(0.5))/(LOG(5x1 0 A 4)-LOG(N'))] 

Fig. 17 Effect of Thermal Fatigue (Cycles) on Thermal Fatigue Strength. 

(Sensitized Model Form) 
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5.7 Model Calibration 

The first estimates of the ultimate and reference values for each effect are given in 
Table 8. First estimates of the empirical material constants, previously determined from linear 
regression of high temperature, mechanical fatigue, creep and thermal fatigue data, are 
summarized in Table 9. These initial estimates were used to calibrate the strength degradation 
model specifically for Inconel 718. Thus, model accuracy is dependent on proper selection of 
ultimate and reference values, which in turn influence the values of the empirical material 
constants. 


Table 8 Initial Estimates for the Ultimate and Reference Values. 


Effect 

Ultimate 

Value 

Symbol 

Estimated 
Ultimate Value 

Reference 

Value 

Symbol 

Estimated 
Reference Value 

Temperature 

Tu 

2369 

To 

75 

Mechanical Fatigue 

Nu 

lxlO 10 

No 

0.5 

Creep 

tu 

1x106 


0.25 

Thermal Fatigue 

N'u 

5x104 

No 

0.5 


Table 9 Initial Estimates for the Empirical Material Constants. 


Effect 

Empirical Material 
Constant Symbol 

Estimated Value 
of Constant 

Applicable 
Temperature (°F) 

High Temperature 

q 

0.2422 

75-1300 

Mechanical Fatigue 

S 

0.3785 

75 

Mechanical Fatigue 

s 

0.2235 

1000 

Mechanical Fatigue 

s 

0.3543 

1200 

Creep 

V 

0.2912 

1000 

Creep 

V 

0.4008 

1100 

Creep 

V 

0.6243 

1200 

Creep 

V 

1.1139 

1300 

Thermal Fatigue 

u 

0.2368 

900 
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As previously mentioned, the quantities used for ultimate and reference values were 
initial estimates. Based on the parameters obtained from linear regression analysis of the data, 
i.e. slope (material constant), y-intercept (log So) and R 2 , an attempt to adjust these initial 
es tima tes to improve the accuracy of the model was made. Noting that the y-intercept value 
corresponds to the log of the reference strength. So, it was necessary to physically define what 
the quantity So represents. For the temperature model, given the data used. So (5.217 or 
164,816 psi) estimates the yield strength of Inconel 7 18 at the reference temperature of 75 °F as 
seen by Figure 3. In order to correlate the So for all effects to the yield strength, the ultimate 
and reference values for mechanical fatigue, creep and thermal fatigue effects were adjusted. 
Adjusting the ultimate value influenced the slope, y-intercept and R 2 values, while adjusting the 
reference value altered the y-intercept value but had no affect on either the slope or R 2 values. 
In addition, certain trends were noted. Increasing the ultimate value increased the So value, 
while increasing the reference value decreased it 

Based on this information, initial estimates were reevaluated for mechanical fatigue, 
creep and thermal fatigue effects. Reevaluation of the initial estimates for the temperature effect 
was not necessary since this temperature data consisted of yield strength values at various 
temperatures, thus So is already correlated to a yield strength value of Inconel 718. For the 
mechanical fatigue effect. Figure 6 shows log So values of 5.1974 (157,543 psi), 5.1067 
(127,850 psi) and 5.1184 (131,341 psi) for temperatures of 75, 1000 and 1200 °F, 
respectively. According to average yield strength data for Inconel 718 [15], these values are 
too low. Therefore, in order to increase these y-intercept values, the ultimate value was varied 
between lxlO 10 and lxlO 11 cycles, while the reference value was varied between 0.5 and 
0.25 cycles. The result was that an ultimate value of lxlO 10 combined with a reference value 
of 0.25 yielded y-intercept values closest to the average yield strength for corresponding 
temperatures. Initial ultimate and reference values for the creep and thermal fatigue models 
were also adjusted accordingly. Figures 18, 19 and 20, show the improved ultimate and 
reference values selected and display the subsequent new linear regression results of the 
mechanical fatigue, creep and thermal fatigue data, respectively. Table 10 lists the improved 
estimates obtained for the ultimate and reference values, while Table 11 provides the 
corresponding new empirical material constants. 
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Figure 1 8 Effect of Mechanical Fatigue (Cycles) cm Fatigue Strength for Inconel 718. 
(Sensitized Model Form Using Improved Estimates) 



Figure 19 Effect of Creep Time (Hours) on Rupture Strength for Inconel 718. 
(Sensitized Model Form Using Improved Estimates) 
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LOG[(LOG(5x10 A 4)- LOG (0.25 ))/( LOG (5x10 A 4)- LOG (N‘))] 

Figure 20 Effect of Thermal Fatigue (Cycles) on Thermal Fatigue Strength. 
(Sensitized Model Form Using Improved Estimates) 
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Table 10 improved Estimates for the Ultimate and Reference Values. 


Effect 

Ultimate 

Value 

Symbol 

Estimated 
Ultimate Value 

Reference 

Value 

Symbol 

Estimated 
Reference Value 

Temperature 

Tu 

2369 

To 

75 

Mechanical Fatigue 

Nu 

lxlO 10 

No 

0.25 

Creep 

tu 

lxlO 5 

to 

0.25 

Thermal Fatigue 

N'u 

5x10* 

N'o 

0.25 


Table 1 1 Improved Estimates for the Empirical Material Constants. 


Effect 

Empirical Material 
Constant Symbol 

Estimated Value 
of Constant 

Applicable 
Temperature (°F) 

High Temperature 

q 

0.2422 

75-1300 

Mechanical Fatigue 

S 

0.3785 

75 

Mechanical Fatigue 

S 

0.2235 

1000 

Mechanical Fatigue 

s 

0.3543 

1200 

Creep 

V 

0.1737 

1000 

Creep 

V 

0.2245 

1100 

Creep 

V 

0.4136 

1200 

Creep 

V 

0.7556 

1300 


Thermal Fatigue 


u 


0.1908 


900 




















CHAPTER 6 

ESTIMATION OF EMPIRICAL MATERIAL 
CONSTANT VARIABILITY 


Due to a lack of sufficient data from which to evaluate the material constants, aj, 
methodology to estimate the variability of these constants was developed. This methodology 
yields estimates for the standard deviations of the constants. For instance, when modeling 
high temperature effects, the material strength degradation model for Inconel 718 is given 
below by equation (6a). 



or 


S = S G 



(14a) 


Taking the log of both sides yields equation (14b) below. 


Log S = -q 


Log 


Tu~T c 

Tu-T 


1 \ 


+ Log S G 


(14b) 


jy 


It is clearly seen that equation (14b) is a linear equation with slope, -q, and y-intercept, 
Log So- Using the temperature data presented in Chapter 5, the linear relationship given by 
equation (14b) is shown graphically in Figure 21. 

Linear regression of this temperature data yielded two parameters, the slope (-0.2422) 
and the y-intercept (5.2170). As previously discussed, the slope was used as a first estimate of 
the empirical material constant for the temperature degradation model. Due to limited 
temperature data, only five data points, concern over the accuracy of this es timated value was 
warranted. Therefore, steps were taken to model this material constant as a random variable so 
that an estimate of its standard deviation could be calculated. 
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Figure 21 Linear Regression of Temperature Data. 


First, maximum and mi nimum feasible slopes and y-intercepts were determined from 
consideration of the data and the linear regression results, such that these extreme parameters 
would theoretically enclose or envelope all actual data. Figure 22 shows the linear regression 
of the temperature data along with postulated maximum and minimum slopes. These extreme 
parameters were obtained by adjusting the slope of the linear regression fit Rotating about the 
y-intercept value, the regression line was adjusted to pass through the outer most points, 
resulting in maximum and minimum slopes. Figure 23 shows the linear regression of the 
temperature data along with maxi mum and minimum y-intercepts. These extreme parameters 
were obtained by shifting the regression line vertically. While maintaining the slope, the 
regression line was shifted to pass through the outer most points, resulting in maximum and 
minim um y-intercept values. 



LOG YIELD STRENGTH (RSI) LOG VELD STRENGTH (PSI) 
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Figure 22 Postulated Maximum and Minimum Slopes. 



0.2421 5x 


0.24215x 


I 


Figure 23 Postulated Maximum and Minimum Y-intercepts. 
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Using the values of the parameters obtained from linear regression along with the 
extreme maximum and min imum values, random variables for slope (-q) and y-intercept 
(log So) were constructed. These random parameters or variables were assumed to have 
normal distributions, with mean values given by the linear regression fits in Figure 21. 
Standard deviation values for the slope and y-intercept were determined using the extreme 
values together with the empirical rule. According to this rule, for a normal distribution, the 
mean value (p) plus or minus three standard deviations (±3o) will contain 99.73% of the 
values [17, 20]. Therefore, the range of the values (maximum value minus the minimum 
value) divided by six yields the standard deviation, a. Although the mean value resulting from 
linear regression (Figure 21) is not equal to p (one-half the range) due to the nature of the data 
and the extreme values obtained, this method provides for an approximation of the standard 
deviation. 



Figure 24 Probability Density Function of a Normal Distribution. 


Values for the standard deviation of the random parameters, slope and y-intercept, 
were estimated as follows: 


'slope 


maximum slope - minimum slope _ 0.2614 — 0.2085 _ ^ ^ 


0.0088 


max imum y - int. - minimum y-int. 167,707.20-162,416.67 001 __ 

Oy-int = = = 88175 


(psi) 


6 


6 
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These random parameters, now expressed in terms of their mean and standard deviation, were 
used to define the probabilistic material strength degradation model for temperature as a random 
parameter model having the following form: 



where -q and So are now random variables for the slope and y-intercept, respectively. 

In order to demonstrate this methodology, modifications were made to PROMISS 
[6]. These modifications included providing random variable input mechanisms for So in 
terms of its mean and standard deviation, adding random number generation capability for So, 
and providing coding to calculate equation (14c), so that results are given in terms of strength, 
S, rather than lifetime strength, S/So- The resulting values for S were calculated by simulation 
using an augmented version of PROMISS called CALLIE92T. Forty values of strength, S, 
corresponding to each temperature value, T, were obtained. Figure 25 displays selected 
strength values of the forty calculated, along with the actual temperature data and the postulated 
envelope of the random parameter model as defined by the extreme parameter values. The 
statistical frequency with which calculated values of S fell within the envelope were noted. 
Since an overwhelmingly large number of S values were found to lie within the envelope, it 
was ascertained that experimental temperature data beyond the known five data points would 
also fall within the envelope. Thus, this estimated value of the standard deviation, rather than 
expat opinion or an assumed value, can be used with greater confidence in the probabilistic 
material strength degradation model embodied in PROMISS. 



Figure 25 Postulated Envelope of Actual and Simulated Temperature (°F) Data. 



CHAPTER 7 

PROBABILISTIC LIFETIME STRENGTH SENSITIVITY STUDY FOR 
MECHANICAL FATIGUE, CREEP AND THERMAL FATIGUE 


A modified version of PROMISS, entitled PROMISS93, was developed for 
sensitizing the model for mechanical fatigue, creep and thermal fatigue effects. Using the 
sensitized probabilistic material strength degradation model embodied in PROMISS93, a 
lifetime strength sensitivity study was conducted. Three effects were included in this study, 
mechanical fatigue, creep and thermal fatigue. The temperature effect was not explicitly 
included as a fourth effect since the data used in this study for the other effects resulted from 
tests conducted at elevated temperatures of 900 to 1000 °F. Therefore, the effect of high 
temperature is inherent in the mechanical fatigue, creep and thermal fatigue empirical material 
constants used to calibrate the model. 

The general form of the multifactor equation given by equation (1), when modified 
for combined mechanical fatigue, creep and thermal fatigue effects, becomes. 
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(15b) 


By making the necessary log transformations to increase model sensitivity and accuracy for 
these three specific effects, equation (15b) becomes. 


s 

log(Nu) - log(No) ~ 

Tlog(tu)-log(to) 

-V 

log(Nu) — log(No) 

So 

log(Nu)-log(N) _ 

_ log(tjj ) log(t) _ 

log(Nu ) — log(N ) _ 


(16a) 


Substitution of the improved ultimate and reference estimates results in equation (16b) below. 


A- 

log(l0 10 ) - log(0.25)" 

-s 

log(l0 5 ) - log(0.25)" 

-V 

log{5xl0 4 )-log(0.25)* 

So 

log{l0 10 )-log(N) 


log(l0 5 )-log(t) 


log(5xl0 4 )-log(N) 


(16b) 
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The ultimate and reference values in equation (16b) became model parameters or 
constraints for the multifactor equation when modified for Inconel 718. Figure 26 illustrates 
these model parameters graphically, wherein each axis represents an effect 

THERMAL FATIGUE (CYCLES) 


MECHANICAL FATIGUE (CYCLES) 

CREEP (HOURS) 

Fig. 26 Inconel 718 Model Parameters for Mechanical Fatigue, 

Creep and Thermal Fatigue Effects. 

Typical sets of input values for the PROMISS model represented by equation (16b) 
are given in Tables 12, 13 and 14. For example. Table 12 shows PROMISS input data for a 
temperature of 1000 °F, a current value of 2.5X10 5 mechanical fatigue cycles, a current value of 
1000 creep hours, and a current value of 2000 thermal fatigue cycles. As seen in Tables 12 
through 14, the above-mentioned current values remain the same with the exception of the 
current value of mechanical fatigue cycles, N. In Tables 13 and 14 the current value of 
mechanical fatigue cycles has been increased to l.OxlO 6 and 1.75X10 6 , respectively. By 
holding two of the three sets of current values constant, sensitivity of lifetime strength towards 
the third set of values, in this case mechanical fatigue cycles, can be ascertained. The complete 
set of current values that were used as input data for this sensitivity study are given in 
Table 15. Notice that the first three rows of the table correspond to the current values listed in 
Tables 12, 13 and 14, respectively. The next three rows of Table 15 show how the current 
values of creep hours were varied, while the last three rows show how the current values of 
thermal fatigue cycles were varied. The results of this study, in the form of cumulative 
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distribution functions, are given in Figures 27 through 29. Figure 27 shows the effect of 
mechanical fatigue cycles on lifetime strength, while Figures 28 and 29 show the effect of 
creep hours and thermal fatigue cycles on lifetime strength, respectively. Note that the c.d.f. 
shifts to the left, indicating a lowering of lifetime strength, as mechanical fatigue cycles 
increase. In this manner, results, in the form of c.d.f.'s, display the sensitivity of lifetime 
strength to any current value of an effect. 

Table 12 Sensitivity Study Input to PROMISS93 for Inconel 718; 

Temperature = 1000 °F and N=2.5xl0 5 Cycles. 


Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(V alue), (% of Mean) 

Mechanical 

Nu 

cycles 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycles 

Normal 

2.5x10 s 

2.5X10 4 

10.0 


No 

cycles 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.2235 

0.0067 

3.0 

Creep 

tu 

hours 

Normal 

1.0x10 s 

l.OxlO 4 

10.0 


t 

hours 

Normal 

l.OxlO 3 

l.OxlO 2 

10.0 


to 

hours 

Normal 

0.25 

0.025 

10.0 


V 

dimensionless 

Normal 

0.1737 

0.0052 

3.0 

Thermal 

N’u 

cycles 

Normal 

5.0X10 4 

5.0X10 3 

10.0 

Fatigue 

N' 

cycles 

Normal 

2.0X10 3 

2-OxlO 2 

10.0 


No 

cycles 

Normal 

0.25 

0.025 

10.0 


u 

dimensionless 

Normal 

0.191 

0.0057 

3.0 


Table 13 Sensitivity Study Input to PROMISS93 for Inconel 718; 
Temperature = 1000 °F and N= l.OxlO 6 Cycles. 

Effect 

Variable Units 

Symbol 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycles 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycles 

Normal 

l.OxlO 6 

1.0x10 s 

10.0 


No 

cycles 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.2235 

0.0067 

3.0 

Creep 

tu 

hours 

Normal 

1.0x10 s 

l.OxlO 4 

10.0 


t 

hours 

Normal 

l.OxlO 3 

l.OxlO 2 

10.0 


to 

hours 

Normal 

0.25 

0.025 

10.0 


V 

dimensionless 

Normal 

0.1737 

0.0052 

3.0 

Thermal 

N'u 

cycles 

Normal 

5-OxlO 4 

5.0x103 

10.0 

Fatigue 

N' 

cycles 

Normal 

2.0x103 

2.0X10 2 

10.0 


N'o 

cycles 

Normal 

0.25 

0.025 

10.0 


u 

dimensionless 

Normal 

0.191 

0.0057 

3.0 







Table 14 Sensitivity Study Input to PROMISS93 for Inconel 718; 
Temperature = 1000 T and N=1.75xl0 6 Cycles 


Effect 

Variable Units 

Symbol 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycles 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycles 

Normal 

1.75x106 

1.75x10 s 

10.0 


No 

cycles 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.2235 

0.0067 

3.0 

Creep 

tu 

hours 

Normal 

1.0x10 s 

l.OxlO 4 

10.0 


t 

hours 

Normal 

l.OxlO 3 

l.OxlO 2 

10.0 


to 

hours 

Normal 

0.25 

0.025 

10.0 


V 

dimensionless 

Normal 

0.1737 

0.0052 

3.0 

Thermal 

N'u 

cycles 

Normal 

5-OxlO 4 

5.0X10 3 

10.0 

Fatigue 

N' 

cycles 

Normal 

2.0xlCP 

2.0X10 2 

10.0 


No 

cycles 

Normal 

0.25 

0.025 

10.0 

——-saas—a— 

u 

dimensionless 

Normal 

0.191 

0.0057 

3.0 


Table 15 Selected Current Values for Sensitivity Study of tbe Probabilistic 
Material Strength Degradation Model for Inconel 718. 


Mechanical Fatigue 
(Cycles) 

Creep 

(Hours) 

Thermal Fatigue 
(Cycles) 

2.5 x 10 5 " 


itoo " 

1.0 x 10« 

1000 

2000 

1.75 x 10* 

1000 

2000 

1.0 x 10 6 

250 

2000 

1.0 x 10 6 

1000 

2000 

1.0 x 10 6 

1750 

2000 

1.0 x 106 

1000 

500 

1.0 x 106 

1000 

2000 

1.0 x IQ 6 

1000 

3500 
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0.40 0.45 0.50 0.55 0.60 0.65 

LIFETIME STRENGTH, S/So 


Fig. 27 Comparison of Various Levels of Uncertainty of Mechanical Fatigue (Cycles) 
on Probable Strength for Inconel 718 for 2000 Thermal Fatigue Cycles 
and 1000 Hours of Creep at 1000 °F. 



0.40 0.45 0.50 0.55 0.60 0.65 

LIFETIME STRENGTH, S/So 

Fig. 28 Comparison of Various Levels of Uncertainty of Creep Time (Hours) 
on Probable Strength for Inconel 718 for lxlO 6 Mechanical Fatigue Cycles 
and 2000 Thermal Fatigue Cycles at 1000 °F. 
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LIFETIME STRENGTH, S/So 

Fig. 29 Comparison of Various Levels of Uncertainty of Thermal Fatigue (Cycles) 
on Probable Strength for Inconel 718 for lxlO 6 Mechanical Fatigue Cycles 
and 1000 Hours of Creep at 1000 °F. 



CHAPTER 8 

MODEL VERIFICATION STUDY 


Using the probabilistic material strength degradation model embodied in PROMISS, a 
model verification study was conducted. The basic assumption, that two or more effects acting 
on the material multiply (i.e., independent variables), was evaluated. Available data allowed 
for a verification study comparing a combination of mechanical fatigue effects at 75 °F and 
temperature effects at 1000 °F to mechanical fatigue effects at 1000 °F. That is, a combination 
of mechanical fatigue and temperature by model was compared to the combination of these two 
effects by experiment. The input values for the combination of these two effects by model are 
given in Tables 16 through 18, while the input values for the combination of these two effects 
by experiment are provided in Tables 19 through 21. Three different current values of 
mechanical fatigue cycles were used so that the verification study would encompass a range of 
fatigue cycle values. The results of this study, in the form of cumulative distribution functions, 
are given in Figures 30 through 32. Figure 30 displays lifetime strength predictions for the 
combination of mechanical fatigue and temperature by model, while Figure 31 displays results 
for the combination of these two effects by experiment. Figure 32 is an overlay of the two sets 
of results. It is evident that there is approximately a 20% difference between the two sets of 
distributions. 

Due to the questionable mechanical fatigue material constant (s = 0.37848) used in the 
combination by model input, a second verification study was conducted. Once again, a 
combination of these two effects by model was compared to the combination by experiment. 
However, an adjusted mechanical fatigue material constant (s = 0.141) was input in place of 
the questionable mechanical fatigue material constant at a temperature of 75 °F. This value was 
estimated by noting the percent difference (37 %) between the calculated slopes at 1000 °F and 
1200 °F. The improved input values for this second verification study are provided in 
Tables 22 through 24. The input values for combination by experiment were the same as 
before. The results are given by Figures 33 through 36. Figure 33, overlays the results for the 
combination by model and those by experiment The 20% difference was greatly reduced. For 
clarity. Figures 34, 35 and 36 overlay the results for both model and experiment for current 
mechanical fatigue cycle values of 2.5x10 s , lxlO 6 and 1.75xl0 6 cycles, respectively. A 
percent difference of less than 5% was observed for all three current mechanical fatigue cycle 
values. 
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Table 16 Verification Study Input to PROMISS93 for Inconel 718; 
Combination by Model, N=2.5xlO s cycles. 
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Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycle 

Normal 

2.5x10 s 

2.5X10 4 

10.0 

(at 75 °F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.3785 

0.0114 

3.0 

High 

Tu 

Y 

Normal 

2369.0 

236.90 

10.0 

Temperature 

T 

Y 

Normal 

1000.0 

100.00 

10.0 

(at 1000 °F) 

To 

Y 

Normal 

75.0 

7.50 

10.0 


q 

dimensionless 

Normal 

0.2422 

0.0088 

3.6 


Table 17 Verification Study Input to PROMISS93 for Inconel 718; 
Combination by Model, N=1.0xl0 6 cycles. 

Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycle 

Normal 

l.OxlO 6 

1.0x10 s 

10.0 

(at 75 °F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s dimensionless 

Normal 

0.3785 

0.0114 


High 

Tu 

Y 

Normal 

2369.0 

236.90 

10.0 

Temperature 

T 

Y 

Normal 

1000.0 

100.00 

10.0 

(at 1000 °F) 

To 

Y 

Normal 

75.0 

7.50 

10.0 


q dimensionless 

Normal 

0.2422 

0.0088 

3.6 


Table 18 Verification Study Input to PROMISS93 for Inconel 718; 


Combination by Model, N=1.75xl0 6 cycles. 

Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycle 

Normal 

1.75X10 6 

1.75x10 s 

10.0 

(at75°F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s dimensionless 

Normal 

0.3785 

0.0189 

3.0 

High 

Tu 

Y 

Normal 

2369.0 

236.90 

10.0 

Temperature 

T 

Y 

Normal 

1000.0 

100.00 

10.0 

(at 1000 °F) 

To 

Y 

Normal 

75.0 

7.50 

10.0 

_ 

q dimensionless 

Normal 

0.2422 

0.0088 

3.6 


i 








Table 19 Verification Study Input to PROMISS93 for Inconel 718; 
Combination by Experiment, N=2.5xlO s cycles. 
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Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

1.0x1010 

1.0x109 

10.0 

Fatigue 

N 

cycle 

Normal 

2.5x10 s 

2.5x10* 

10.0 

(at 1000 °F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.2235 

0.0067 



Table 20 Verification Study Input to PROMISS93 for Inconel 718; 
Combination by Experiment, N=1.0xl0 6 cycles. 


Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

1.0x109 

10.0 

Fatigue 

N 

cycle 

Normal 

l.OxlO 6 

1.0x10 s 

10.0 

(at 1000 °F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.2235 

0.0067 

3.0 


Table 21 Verification Study Input to PROMISS93 for Inconel 718; 
Combination by Experiment, N=1.75xl0 6 cycles. 


Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycle 

Normal 

1.75x106 

1.75x10 s 

10.0 

(at 1000 °F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.2235 

0.0067 

3.0 
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0.50 0.55 0.60 0.65 0.70 0.75 

LIFETIME STRENGTH, S/So 

Figure 30 Comparison of Various Levels of Uncertainty of Mechanical Fatigue (Cycles) 

on Probable Strength for Inconel 718. 

(Combination of Mechanical Fatigue and High Temperature Effects by Model) 



0.70 0.75 0.80 0.85 0.90 0.95 

LIFETIME STRENGTH, S/So 

Figure 31 Comparison of Various Levels of Uncertainty of Mechanical Fatigue (Cycles) 

on Probable Strength for Inconel 718. 

(Combination of Mechanical Fatigue and High Temperature Effects by Experiment) 
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LIFETIME STRENGTH, S/So 


Figure 32 Overlay of Results for the Combination of Mechanical Fatigue and Temperature 

Effects by Model and Experiment 
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Table 22 Modified Verification Study Input to PROMISS for Inconel 718; 
Combination by Model, N =2.5x105 cycles. 


Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

l.OxlO^ 

10.0 

Fatigue 

N 

cycle 

Normal 

2.5x105 

2.5X10 4 

10.0 

(at 75 °F) 

N 0 

cycle 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.141 

0.0042 

3.0 

High 

Tu 

Y 

Normal 

2369.0 

236.90 

10.0 

Temperature 

T 

Y 

Normal 

1000.0 

100.00 

10.0 

(at 1000 °F) 

To 

Y 

Normal 

75.0 

7.50 

10.0 


q 

dimensionless 

Normal 

0.2422 

0.0088 

3.6 


Table 23 Modified Verification Study Input to PROMISS for Inconel 718; 
Combination by Model, N=1.0xl0 6 cycles. 


Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(Value), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycle 

Normal 

l.OxlO 6 

l.OxlO 5 

10.0 

(at75°F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.141 

0.0042 

3.0 

High 

Tu 

Y 

Normal 

2369.0 

236.90 

10.0 

Temperature 

T 

°F 

Normal 

1000.0 

100.00 

10.0 

(at 1000 °F) 

To 

Y 

Normal 

75.0 

7.50 

10.0 


q 

dimensionless 

Normal 

0.2422 

0.0088 

3.6 


Table 24 Modified Verification Study Input to PROMISS for Inconel 718; 
Combination by Model, N=1.75xl0 6 cycles. 


Effect 

Variable 

Symbol 

Units 

Distribution 

Type 

Mean 

Standard Deviation 
(V alue), (% of Mean) 

Mechanical 

Nu 

cycle 

Normal 

l.OxlO 10 

l.OxlO 9 

10.0 

Fatigue 

N 

cycle 

Normal 

1.75X10 6 

1.75x105 

10.0 

(at 75 °F) 

No 

cycle 

Normal 

0.25 

0.025 

10.0 


s 

dimensionless 

Normal 

0.141 

0.0042 

3.0 

High 

Tu 

Y 

Normal 

2369.0 

236.90 

10.0 

Temperature 

T 

Y 

Normal 

1000.0 

100.00 

10.0 

(at 1000 °F) 

To 

Y 

Normal 

75.0 

7.50 

10.0 


q 

dimensionless 

Normal 

0.2422 

0.0088 

3.6 
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LIFETIME STRENGTH, S/So 


Figure 33 Overlay of Results for the Combination of Mechanical Fatigue and Temperature 
Effects by Model (Using Estimated Value of s) and Experiment 



LIFETIME STRENGTH, S/So 


Figure 34 Overlay of Results for the Combination of Mechanical Fatigue and Temperature 
Effects by Model (Using Estimated Value of s) and Experiment; N=2.5xl0 5 Cycles. 
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LIFETIME STRENGTH, S/So 

Figure 35 Overlay of Results for the Combination of Mechanical Fatigue and Temperature 
Effects by Model (Using Estimated Value of s) and Experiment; N=1.0xl0 6 Cycles. 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

0 

LIFETIME STRENGTH, S/So 

Figure 36 Overlay of Results for the Combination of Mechanical Fatigue and Temperature 
Effects by Model (Using Estimated Value of s) and Experiment; N=1.75xl0 6 Cycles. 




CHAPTER 9 
DISCUSSION 


To ensure model accuracy in lifetime strength predictions, close attention was paid to 
model sensitization and calibration. When the current value and the reference value were small 
compared to the ul tima te value, model transformation, by taking the log of each value within 
the product term, was required for model sensitivity. As shown for mechanical fatigue, creep 
and thermal fatigue effects in Figures 5 through 6, 8 through 9, and 16 through 17, 
respectively, this transformation resulted in considerable increases in the linear regression R 2 
values. The closer the R 2 value is to a value of one, the better the linear regression fit 

Calibration of the model specifically for Inconel 718 required actual experimental 
data. Based on this data, initial ultimate and reference values for each effect were estimated and 
are provided in Table 8. Linear regression of data individually for each effect resulted in initial 
estimates for the empirical material constants. These constants for temperature, mechanical 
fatigue, creep and thermal fatigue effects are given in Table 9. Further calibration involved 
adjusting these initial estimates so that y-intercept (log So) values, resulting from linear 
regression analysis, corresponded to average yield strength values of Inconel 718 at specified 
temperatures. By correlating the So values for all effects to average yield strengths, accuracy 
in modeling two or more effects was increased. These improved estimates are given in 
Tables 10 and 1 1. These estimates were used for the mean values in the sensitivity study input 
files (Tables 12 through 14) to PROMISS93. 

Methodology for estimating the variability of the empirical material constants was 
developed in Chapter 6 as a means for dealing with limited data. For the temperature effect, a 
standard deviation value of 0.0088 or 3.6% of the mean slope (0.2422) was calculated. This 
value, rather than expert opinion, may be used with greater confidence in the probabilistic 
material strength degradation model embodied in PROMISS93. Parallel steps may be taken to 
determine standard deviation estimates for the empirical material constants of the other effects. 

The sensitivity study, discussed in Chapter 7.0, included only three effects, 
mechanical fatigue, creep and thermal fatigue, as modeled by equation (16b). The results of 
this study, in the form of cumulative distribution functions, are given in Figures 27 through 
29. The sensitivity of lifetime strength to the number of mechanical fatigue cycles is seen by 
the shift of the c.d.f. to the left in Figure 27 as the number of cycles increases from 2.5X10 5 to 
1.75X10 6 . The same phenomenon is seen in Figures 28 and 29. Thus, increasing the current 
number of the variable decreased the predicted lifetime strength as expected. The temperature 
effect was not explicitly included in this study due to the fact that data for the other three effects 
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resulted from tests conducted in a high temperature environment (900 °F to 1000 °F). Thus, 
the effect of temperature is inherent in the estimated empirical material constants for the other 
three effects. This is evidenced by the changing slopes in Figure 19 for the creep effect The 
slope or material constant changes according to the test temperature. At a test temperature of 
1000 °F, the material constant (slope) is -0.17372, but increases with temperature to a "steeper" 
value of -0.75557 at a test temperature of 1300 °F. An increase in the material constant with an 
increase in temperature is expected. However, as seen by Figure 18, the mechanical fatigue 
material constant (slope) is highest at the lowest test temperature of 75 °F. Since this slope is 
based upon only four questionable data points, it is presumed to be inaccurate. Therefore, 
based on observed trends in the change of slopes for the mechanical fatigue effect at 
temperatures of 1000 °F and 1200 °F (Figure 18), an adjusted value for the mechanical fatigue 
material constant at 75 °F was determined. The result was a modified slope 37% less than the 
slope obtained at a temperature of 1000 °F. Without additional mechanical fatigue data at a test 
temperature of 75 °F, this adjusted slope can be neither confirmed nor rejected. 

Both the questionable (s = 0.37848) and the adjusted (s = 0.141) mechanical fatigue 
material constants at 75 °F were used in verification studies presented in Chapter 8. Available 
data allowed for a verification study comparing a combination of mechanical fatigue and 
temperature effects by model to the combination of these two effects by experiment. The 
results of this study, in the form of c.tLf.'s, are given in Figures 30 through 32. The 
sensitivity of lifetime strength to the number of current mechanical fatigue cycles is seen by the 
shift of the c.d.f. to the left (Figures 30 and 31) as the number of cycles increases. Thus, 
increasing the number of current fatigue cycles decreases the predicted lifetime strength as 
expected. As seen by the overlay of distributions in Figure 32, there is approximately a 20% 
difference between the results obtained by model and those obtained by experiment A major 
possibility for this large discrepancy is the questionable mechanical fatigue material constant at 
75 °F. To test this assumption, a second parallel verification study using the adjusted 
mechanical fatigue material constant value was conducted. The results are given in Figures 33 
through 36. Comparison of Figure 36 to Figure 33 shows a substantial decrease in the 
discrepancy between the two sets of distributions. From Figures 34 through 36, the percent 
difference between the results is less than 5% for all three current values of fatigue cycles 
evaluated. Thus, the questionable mechanical fatigue material constant calculated from the 
mechanical fatigue data at 75 °F was responsible for a large percent of the discrepancy between 
the initial results from the first verification study. 
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CHAPTER 10 
CONCLUSIONS 


A probabilistic material strength degradation model, applicable to aerospace materials, 
has been postulated for predicting the random lifetime strength of structural components for 
propulsion system components subjected to a number of effects. This model, in the form of a 
randomized multifactor equation, has been developed for four effects, namely, high 
temperature, mechanical fatigue, creep and thermal fatigue. Inconel 718 data for these effects 
was obtained from the open literature. Based on this data, initial ultimate and reference values 
were estimated. It was determined that when the current and reference values are small 
compared to the ultimate value the model is insensitive. Therefore, a transformation to 
sensitize the model for the effects of mechanical fatigue, creep and thermal fatigue was 
required. Model transformation resulted in significant increases in the R 2 (goodness of fit) 
values. The current version of PROMISS, entitled PROMISS93, provides for this 
transformation for these three effects. 

Linear regression of the data for each effect resulted in estimates for the empirical 
material constants, as given by the slope of the linear fit. These estimates, together with 
ultimate and reference values, were used to calibrate die model specifically for Inconel 718. By 
adjusting these initial estimates so that the y-intercept or So values corresponded to average 
yield strength values of Inconel 718, accuracy in modeling two or more effects was improved. 
Thus, model accuracy is dependent on the proper selection of ultimate and reference values, 
which in turn influence the values of the empirical material constants used in calibration of the 
model. Calibration of the model for other materials is also dependent on experimental data and 
is not possible without it 

Methodology for estimating the standard deviation of empirical material constants 
offered a way for dealing with limited data. This methodology results in better estimates of the 
standard deviations based on actual experimental data, rather than expert opinion. Lack of 
sufficient data from which to evaluate the material constants warranted the development of this 
methodology. 

Results from a sensitivity study involving mechanical fatigue, creep and thermal 
fatigue effects showed that the c.d.f.'s shift to the left, indicating a lowering of lifetime 
strength, for increasing current values of an effect Further development and evaluation of this 
three effect model, as well as other models, requires that it be compared to real responses of 
Inconel 718 samples subjected to these combined effects during experimentation. Thus, 
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additional experimental data is crucial for the continued development and evaluation of the 
probabilistic material strength degradation model presented in this thesis. 

Limited verification studies involving two effects, mechanical fatigue arid high 
temperature, were conducted. Results showed a combination of the two effects by model to be 
more conservative than the combination by experiment The first verification study yielded a 
20% discrepancy between the results obtained by model and those obtained by experiment. 
Questionable mechanical fatigue data at a temperature of 75 °F is presumed to be a major cause 
of the discrepancy. This conclusion was drawn after conducting a second verification study 
using an adjusted value in place of the questionable one. The outcome was a significant 
reduction in the discrepancy, from 20% to less than 5%, between the results of a combination 
of these two effects by model and the combination by experiment Therefore, the data, rather 
than the nature of the model, is the presumed source of error. Thus, the basic assumption of 
the model, that two or more effects multiply (i.e., effects are independent), is strongly 
supported by this limited verification study. The remaining 5% difference may be due to the 
lack of uniformity among the specimens tested. As seen by Table A.5 in the Appendix, 
specimen shape and heat treatment varied between the effects. Specimen shape, as well as heat 
treatment, can influence material properties. Another reason for the 5% difference may be 
synergistic effects (i.e., dependence between effects). As previously discussed, equation (1) is 
an approximated solution to a separable partial differential equation. In order to account for 
synergistic effects and perhaps eliminate this 5% difference, additional terms would have to be 
added to equation (1 ). The resulting reduction in error may or may not warrant complication of 
the model by the inclusion of additional terms. Based on the results obtained from the second 
verification study, this complication is not warranted. However, additional verification studies 
for the combination of other effects must first be conducted before a more refined model can be 
developed. As previously discussed, the availability of experimental data will determine 
whether or not further studies can be conducted. 

In conclusion, methodology for improving lifetime strength prediction capabilities is 
presented. The probabilistic material strength degradation model in the form of a randomized 
multifactor equation is developed for four effects and calibrated to best reflect physical reality 
for Inconel 718. Systematic and repeatable methods of model calibration and evaluation are 
developed. Basic understanding and evaluation of the model is generated through sensitivity 
and verification studies. The sensitivity of random lifetime strength to any current value of an 
effect can be ascertained. Probability statements in the form of cumulative distribution 
functions allow improved judgments to be made regarding the likelihood of lifetime strength, 
thus enabling better design decisions to be made. 
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APPENDIX 


This appendix provides the experimental Inconel 718 data analyzed by the postulated 
material strength degradation model. The purpose of this appendix is to allow the calculations 
of Chapter 5 to be repeated. Data for all effects will be presented in tabular form. Tables A. 1 - 
A.4 present the high temperature, mechanical fatigue, creep and thermal fatigue data, 
respectively. Table A.5 provides reference numbers and figure numbers for displayed data, as 
well as, specimen and heat treatment specification s for all data presented in this thesis. 


Table A. 1 Inconel 718 High Temperature Tensile Data. 


TEST TEMPERATURE, 
Y 

TENSILE STRENGTH, 
PSI 

7.50E+01 

1.63E+05 

6.00E+02 

1.56E+05 

1.00E+03 

1.48E+05 

1.20E+03 

1.40E+05 

1.30E+03 

1.35E+05 


Table A.2 Inconel 718 Mechanical Fatigue Data. 


TEST 

TEMPERATURE, 

Y 

FATIGUE STRENGTH, PSI 

10 5 

CYCLES 

10 6 

CYCLES 

10 7 

CYCLES 

10* 

CYQ.ES 

75 

132,000 

101,000 

92,000 

90,000 

1000 

111,000 

102,000 

95,000 

90,000 

1200 

100,000 

94,000 

88,000 

72,000 
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Table A.3 Inconel 718 Rupture Data. 


TEST 

TEMPERATURE, 

°F 

RUPTURE LIFE, 
HRS 

RUPTURE 

STRENGTH, 

PSI 


1000 

27.8 

158000 



133.2 

150000 



256.0 

145000 



814.9 

140000 



1731.0 

134000 



8473.0 

124000 



21523.6 

118000 


1100 

28.2 

135000 



62.0 

130000 



151.9 

123000 



367.5 

117000 



2327.6 

105000 



10606.2 

94000 



33990.7 

86000 


1200 

10.6 

115000 



30.8 

108000 



150.0 

96000 



747.2 

87000 



3131.5 

78000 



7263.0 

68000 



10232.0 

63000 


1300 

18.0 

86000 



70.5 

76000 



182.7 

68000 



476.8 

60000 



808.0 

55000 



2870.7 

44000 



6048.0 

37000 
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Table A.4 Inconel 718 Thermal Fatigue Data. 


Cycles to Failure, 
N'f 

Reversals to 
Failure, 2NV 

Total Strain 
Amplitude 

Plastic Strain 
Amplitude 

45 

90 

0.01 

0.005 

140 

280 

0.0075 

0.0029 

750 

1500 

0.005 

0.0011 

9750 

19500 

0.004 

0.0003 


Table A.5 Inconel 718 Data Summary. 


EFFECT 

REFERENCE 

NUMBER 

FIGURE 

NUMBER 

SPECIMEN 

HEAT TREATMENT 

Temperature 

[14] 

2, 3, 21, 
22, 23, 25 

hot-rolled round, 
4-inch diameter, 
from single sheet 

1950°F/1 hr, plus 
1400°F/10 hr, F.C. 

100 °F/hr to 1200°F, 
hold at 1200°F for 8 hr 

Mechanical 

Fatigue 

[14] 

4, 5, 6, 18 

forging, 
hot-rolled bar, 
average grain 
size of 0.0008 in 

1750°F/1 hr, plus 
1325°F/8 hr, F.C. to 
1 1 50°F, hold at 1 1 50°F, 
total aging time of 18 hr 

Creep 

[2] 

7, 8, 9, 19 

flat-pancake, 

21 in diameter x 
1 in thick 

1800°F/2 hr, A.C., plus 
1325°F/8 hr, F.C. 
100°F/hr to 1 150°F/8 hr, 
A.C. 

Thermal 

Fatigue 

[16] 

10, 11, 12, 
13, 14, 15, 
16, 17, 20 

forging, round, 
11 mm diameter, 
gage length of 
15 mm 

1253K x 1 hr, W.Q., 
997K x 8 hr -(55K/hr) 
to 893K x 8 hr, A.C. 
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